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(U)  Predictions  of  the  amount  of  water  thrown  upward  from  shallow  depth 
explosions  are  presented  in  support  of  a  new  water  barrier  concept  for  the 
protection  of  Riverine  vessels.  The  predictions  are  obtained  using  a  model 
based  on  a  generalized  formulation  of  hydrodynamics  using  an  incompressi¬ 
ble  liquid  assumption,  which  has  been  previously  validated  on  a  variety  of 
applications  including  water  plume  predictions.  Depending  on  both  the 
charge  weight  and  initial  depth,  this  study  presents  estimates  for  the  Effective 
Water  Length  (EWL)  values  at  heights  and  time  intervals  relevant  for  riverine 
vessel  protection  from  Rocket  Propelled  Grenade  (RPG)  attacks. 


(U)  INTRODUCTION 

(U)  A  new  water  barrier  concept  is  being  investigated  to  support  riverine  vessel  RPG  survivabil¬ 
ity.  A  water  barrier  could  protect  maritime  vessels  against  a  large  range  of  threats  through  the 
formation  of  an  interposing  wall  of  water.  This  is  accomplished  through  the  delivery  and  detona¬ 
tion  of  an  explosive  charge  beneath  the  waterline  (~2-3ft)  along  the  incident  threat  trajectory  12  ft 
from  a  central  point  above  the  vessel,  forming  a  literal  water  barrier  that  initiates  and  absorbs  an 
incident  shape  charge  jet.  The  water  barrier  system  order  of  operation  includes  an  IR  detection  of 
threat  launch,  radar  trajectory  tracking,  intercept  solution,  engagement,  water  barrier  formation, 
RPG  initiation/jet  formation,  and  jet.  An  additional  benefit  of  this  system  is  the  minimal  collat¬ 
eral  damage  caused  by  the  water  barrier  where  blast  and  fragmentation  are  mitigated  underneath 
the  waterline. 

(U)  The  water  barrier  concept  allows  for  the  use  of  available  maritime  resources  while  relaxing 
the  precision  requirements  for  engagement  timing  and  threat  trajectory  height  determination.  By 
using  a  simplified  2D  threat  trajectory,  the  use  of  low  cost  radar  systems  (less  than  $20,000)  can 
be  enabled.  The  relaxed  engagement  timing  results  from  a  large  duration  of  time  (~  1.0  seconds) 
where  the  water  barrier  remains  of  adequate  size  to  defeat  incident  RPGs.  Initial  shape  charge  jet 
performance  studies  against  water  show  a  requirement  of  48”  water  depth  to  completely  erode  an 
initiated  RPG-7  jet.  Additional  factors  related  to  the  dynamics  during  the  water  barrier  formation 
may  additionally  reduce  the  required  depth  to  erode  the  incident  shape  charge  jet. 
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(U)  Water  barrier  formation  eombined  with  detection/engagement  timing  caleulations  show  that 
the  water  barrier  ean  defend  against  engagements  as  close  as  149  ft  (45  meters)  (shooter  to  plat¬ 
form)  at  a  water  depth  of  2-3  ft.  The  water  barrier  also  removes  fratricide/blast  lethality,  which  is 
a  key  issue  for  lightly  armored  riverine  vessels. 

(U)  It  should  be  noted  that  the  water  barrier  can  be  used  in  all  brown  and  green  water  environ¬ 
ments  as  the  width  of  the  engagement  range  (width  beyond  the  vessel  where  water  barrier  is 
formed  [12  ft])  of  both  sides  of  the  riverine  vessel  is  shorter  than  the  length  of  the  vessel  (The  riv¬ 
erine  vessel  typically  requires  a  draft  of  3  ft  and  river  width  wide  enough  to  turn  the  vessel 
around). 

(U)  Within  the  water  barrier  effort  there  are  currently  multiple  variables  that  must  be  studied  and 
integrated  to  demonstrate  the  effectiveness  in  defeating  an  incident  RPG  threat.  This  paper  will 
focus  on  the  computational  modeling  of  the  water  barrier  formation  through  energetic  initiation. 
Timing  (formation/dwell)  and  Equivalent  Water  Length  (EWL)  will  be  studied  as  functions  of: 
charge  weight,  detonation  depth,  wall  height  and  physical  boundary  conditions.  Additional  studies 
of  the  ability  to  induce  an  upward  angular  velocity  to  the  incident  threat  will  also  be  performed. 


(U)  MODELING  DETAILS 

(U)  The  predictions  shown  in  this  study  are  obtained  using  a  model  based  on  a  generalized  formu¬ 
lation  of  hydrodynamics  using  an  incompressible  liquid  assumption  [1].  This  formulation  is  well 
suited  for  predicting  long-time  bubble  and  plume  dynamics  for  the  several  reasons.  First  of  all,  the 
“water”  or  “liquid”  region  is  modeled  as  incompressible,  thereby  allowing  for  time  steps  propor¬ 
tional  to  the  inverse  of  the  water  velocity  as  opposed  to  the  much  smaller  time  steps  that  a  com¬ 
pressible  formulation  would  require  based  on  the  speed  of  sound  in  water.  This  is  important  be¬ 
cause  plume  behavior  occurs  on  the  order  of  seconds.  Secondly,  the  model  allows  for  regions  of 
“spray,”  which  is  typical  of  plume  behavior  in  which  a  well  defined  interface  between  the  bubble 
and  water  or  especially  the  water  and  the  air  does  not  exist.  Finally,  the  computational  model  uses 
a  fixed  “Eulerian”  grid  providing  for  generality  in  studying  complex  bubble  dynamics  and  free 
surface  topology  changes.  For  shallow-depth  explosions  this  includes  the  underwater  bubble 
forming  one  or  more  annular  regions  as  a  downward  moving  jet  intersects  the  bottom  surface  of 
the  bubble  as  it  collapses,  in  addition  to  the  radial  plumes  ejected  on  the  bubble’s  second  expan¬ 
sion,  and  the  eventual  venting  of  the  bubble  into  the  atmosphere. 

(U)  The  generalized  formulation  of  hydrodynamics  is  based  on  the  solution  of  the  mass  and  mo¬ 
mentum  conservation  equations 


A  +  V*(pu)  =  0 
(pu),  +  V  •  (puu)  =  -pgk  -  Vp 


subject  to  the  constraint 


(1) 

(2) 

(3) 
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where  p  is  the  time  and  spatially  varying  density,  is  the  constant  density  of  the  incompressible 
liquid,  u  is  the  velocity  vector,  p  is  the  pressure,  and  k  is  the  unit  vector  in  the  direction  of 
gravity,  g .  These  equations  are  approximated  using  a  split-step  approach.  First,  the  conservations 
equations  (1)  and  (2)  are  approximated  using  an  explicit  second-order  Godunov  Method  with  mo- 
notonized  slope  limiting  without  regards  to  the  constraint  (3).  The  density  constraint  is  then  im¬ 
posed  through  the  solution  of  a  variational  inequality.  We  remark  that  the  imposition  of  (3)  also 
has  an  effect  on  the  momentum,  otherwise  non-physical  and  numerically  unstable  energy  increas¬ 
es  may  occur.  The  treatment  of  this  constraint  is  one  of  the  unique  features  of  the  generalized  hy¬ 
drodynamics  formulation,  distinguishing  it  from  other  approaches,  which  simply  truncate  the  den¬ 
sity,  e.g.  [2].  Finally,  the  pressure  is  determined  using  a  projection  method  based  on  discretiza¬ 
tions  derived  from  a  finite  element  method. 

(U)  These  algorithms  have  been  implemented  within  the  computation  codes  BUB2D  (for  two  di¬ 
mensional  and  axially  symmetric  problems)  and  BUB3D  for  fully  three-dimensional  problems. 
Validations  of  these  codes  on  a  variety  of  applications  including  water  plume  predictions  may  be 
found  in  [1,3-5],  as  well  as  much  more  detail  of  the  numerical  procedures  used. 


(U)  NUMERICAL  RESULTS 

(U)  Computations  were  performed  in  an  axially  symmetric  domain  with  0  s  r  s  /? ,  and 
Zg  <  z  <  Zj.  where  z  =  0  corresponds  to  the  ambient  air-water  free  surface.  The  grids  used  for  the 
computations  were  comprised  of  uniform  square  cells  of  size  h  in  the  sub-region  0  <  r  <  8 , 
-10<z<10  in  units  of  feet.  Outside  of  the  uniform  region  the  grids  are  stretched  to  the  boundary 
of  the  domain  in  which  /?  =  40,  Zg  =-20  and  Z^.  =120  .  In  order  to  study  the  numerical  errors, 
computations  were  performed  using  a  “coarse”  grid  using  a  total  of  120x400  cells  with  /j  =  0.1 , 
and  a  “fine”  grid  using  240  x  800  cells  with  h  =  0.05  . 

(U)  Effective  Water  Length 

(U)  The  results  are  presented  using  density  contours  at  different  times  of  the  bubble  evolution  as 
well  as  measurements  of  the  equivalent  water  length  as  a  function  of  both  height  above  the  sur¬ 
face  and  time.  For  this  analysis  we  consider  the  idealized  case  of  a  missile  travelling  along  a  hori¬ 
zontal  line  at  a  given  height  z  =  .  Setting  the  initial  location  of  the  charge  at  coordinates 

[x^,y^,z^)  =  {0,0,-d)  and  assuming  that  the  missile  travels  on  a  straight  line  trajectory  from 
{-x„,S,Zh)  moving  in  the  positive  x-direction  with  an  infinite  velocity.  Then  the  Equivalent  Water 
Length  for  this  trajectory  at  time  t  is 


1  °° 

EWL{t,S,ZH) - /  p{t,x,S,z„)dx  (4) 

Po  -Xo 

where  p^is  the  ambient  density  of  the  water,  p{t,x,S,z„)  is  the  density  of  the  spray  (water-air) 
mixture  at  time,  S  is  the  offset  of  the  trajectory  from  the  charge  center,  and  z^  is  the  height  of 
the  trajectory. 
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(U)  The  values  for  p{t,x,S,z„)  will  be  approximated  using  an  axially  symmetrie  eomputation  us¬ 
ing  the  BUB2D  eode  which  produced  values  at  discrete  times  t  =  t",  at  cell  center  locations 

{ri,Zj)  where  r  =  +y^  is  the  distance  from  the  axis.  We  can  also  express 

1  °° 

EWL[t,S,z„)  =  —  f  p,,(t,r(x,S),Zff)dx  (5) 

Po  -x„ 


The  integral  (5)  is  evaluated  under  the  assumption  that  is  piecewise  linear  in  time  and  piece- 
wise  constant  in  space  using  the  computational  values  within  each  grid  cell.  Note  that  the  maxi¬ 
mum  values  of  (5)  will  usually  occur  when  5  =  0  which  will  be  evaluated  as  a  “best  case  scenar¬ 
io”.  The  goal  of  these  predictions  is  to  demonstrate  conditions  under  which  the  £WZ  values  are 
larger  than  4  ft.  at  a  height  of  8  feet,  for  a  sufficiently  large  time  interval. 


(U)  Run  Matrix 


(U)  The  initial  plan  called  for  predictions  using  charge  weights  of  1,  1.5,  and  2  lbs,  at  depths  of  3, 
4  and  5  feet.  An  important  parameter  for  the  behavior  of  shallow  depth  explosion  is  the  scaled 
depth  defined  by 


C  = 


D 


(6) 


The  scaled  depths  corresponding  to  this  set  of  tests  are  listed  in  Table  1.  For  some  critical  value  of 
C  =  C*  <  0.56  the  bubble  gas  will  escape  into  the  atmosphere  during  its  first  expansion  and  con¬ 
traction  phase.  Although  the  critical  value  has  not  been  precisely  defined  in  the  literature  (e.g. 
Kedrinski  [6]  gives  the  value  as  1,  but  Blake  and  Gibson  [7]  clearly  show  an  example  where 
C  =0.56  which  does  not  vent).  It  is  likely  that  the  there  is  a  transition  region  where  the  bubble 
partially  vents  into  the  atmosphere  with  the  amount  of  venting  depending  C.  In  [8],  Young  gives 
the  bound  C*  =0.2  below  which  the  gas  can  be  observed  surrounding  a  central  column  of  water 
in  the  plume.  Experiments  performed  at  the  University  of  Maryland  [9]  using  small-scale  RP80 
charges  on  a  sandy  bottom  yielded  an  approximate  value  C*  =  0.25  in  rough  agreement  with 
Young’s  estimate.  Based  on  that  data  we  expect  little  or  no  venting  even  for  the  range  of  trials 
listed  in  Table  1  where  C  >  0.448 . 


Table  1  —  Scaled  depths  for  the  original  run  matrix. 


Depth  (ft) 

1.0 

Charge  Weight  (lbs) 

1.5 

2.0 

3 

0.565 

0.493 

0.448 

4 

0.760 

0.664 

0.603 

5 

0.959 

0.838 

0.761 

We  remark  that  these  computations  were  extended  to  shallower  depths,  and  these  results  will  ap¬ 
pear  in  a  more  extensive  report,  currently  in  preparation. 

(U)  Grid  effects  and  Numerical  error 
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(U)  We  begin  by  using  the  same  resolution  for  the  grids  described  in  the  previous  section  and  ap¬ 
proximate  the  known  solution  for  the  case  of  a  spherical  bubble  in  an  infinite  incompressible  me¬ 
dium,  governed  by  the  Rayleigh-Plesset  equation  [10]  with  conditions  set  so  that  the  ratio  of  the 
initial  (minmum)  and  maximum  bubble  radius  was  =13.28  and  the  period  of  oscilla¬ 

tion  is  =  0.372236  s. 

(U)  For  this  infinite  domain  problem  there  are  three  sources  of  error,  namely  the  nominal  spatial 
discretization  size  h,  the  time  step  sizes  and  the  boundary  truncation.  The  time  steps  are  selected 
adaptively  based  on  changes  in  the  bubble  volume.  The  values  listed  in  Table  2  compare  the  er¬ 
rors  in  the  bubble  period  (time  to  first  minimum)  and  bubble  maximum  radius  from  two  separate 
computations  of  the  infinite  domain  problem.  In  this  table  the  value  L  is  the  distance  from  the 
computational  boundary  to  the  charge  center,  is  the  number  of  cells  used  in  each  direction  for 
the  computation,  Nj  is  the  total  number  of  time  steps  taken  before  the  final  time  of  =0.4. 
Note  that  the  value  of  on  the  fine  grid  is  more  than  double  the  value  on  the  coarser  grid  since 
the  domain  size  was  also  doubled  in  each  direction.  The  computational  results  are  indicated  by 
for  the  maximum  bubble  radius  and  for  the  first  bubble  period.  The  value  is  deter¬ 
mined  from  the  maximum  volume  computed  from  the  density  field  for  a  sphere  of  equivalent  vol¬ 
ume.  The  errors  (E)  are  listed  as  relative  errors  between  the  computed  values  and  the  exact  values 
listed  above.  These  results  indicate  the  accuracy  of  the  code  for  obtaining  the  critical  bubble  dy¬ 
namics  parameters,  which  in  turn  affect  the  dynamics  of  shallow  water  plumes.  For  the  computa¬ 
tions  of  the  shallow  plume  dynamics  presented  next,  all  the  results  were  based  on  the  /i  =  0.05 
grid  so  that  we  expect  approximately  1%  error  in  the  maximum  equivalent  bubble  radius  and 
about  4%  for  the  period. 


Table  2  —  Symmetric  bubble  simulation  results. 


h 

L 

Ah 

'^max 

rpih 

e(t^) 

0.1 

50 

140 

823 

6.5735 

1.77% 

0.3448 

7.38% 

0.05 

100 

400 

1581 

6.6382 

0.81% 

0.3567 

3.94% 

(U)  Fig.  2  displays  the  difference  from  the  coarse  and  fine  grids  for  the  case  for  a  2  lb.  charge 
weight  and  an  initial  depth  of  3  feet  (W2D3)  at  early  times  t  <  0.5  s  .  This  represents  a  shallow 
depth  non-venting  case.  The  top  row  displaying  the  contours  resulting  from  the  coarse  grid,  and 
the  corresponding  fine  grid  results  immediately  below.  The  results  from  both  grids  are  qualitative¬ 
ly  very  similar,  with  each  predicting  a  central  plume  above  the  bubble  at  t  =  0.1  s  when  the  bubble 
is  near  its  maximum  volume.  As  the  bubble  begins  to  compress  a  downward  moving  plume  is 
seen  piercing  through  the  bubble  and  impacting  the  bottom  surface  at  t  =  0.2  s.  During  the  next 
re-expansion  radial  plumes  are  ejected  surrounding  the  central  plume  as  seen  at  t  =  0.5  s.  The  evo¬ 
lution  of  the  plumes  rising  and  falling  at  later  times  i  a  1  s  are  matched  very  well  on  each  grid  as 
seen  comparing  the  top  (coarse  grid)  and  bottom  (fine  grid)  rows.  The  relative  density  palette 
used  for  the  contours  are  shown  below  the  frame  images  in  both  figures,  so  that  blue  corresponds 
to  pure  water,  and  white  corresponds  to  void,  and  regions  of  “spray”  between  the  two. 

U)  Water  Barrier  Measurements 

(U)  For  the  water  barrier  measurements  we  computed  the  values  EWL{t,0,z)  given  by  (5) 

through  the  center  of  the  plume  using  a  density  contour  plot  shown  in  Fig.  3  for  the  D3W2  case 
using  the  fine  grid.  Below  the  EWL  contour  are  the  corresponding  density  contours  at  various  dis- 
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Crete  times  and  heights  from  0  to  10  feet  above  the  surface.  In  the  contour  plot,  EWL  values  larger 
than  5  are  simply  shown  as  light  red  as  indicated  in  the  palette  on  the  right. 


0.0  0.2  0.4  0.6  0.8  1.0 

Fig.  2.(U)  Run  W2D3  coarse  (top)  and  fine  (bottom)  grid  density  contours  at  early  times. 


(U)  The  contours  of  EWL  values  for  the  original  run  matrix  in  Table  1  are  shown  in  Fig.  4.  This 
figure  clearly  indicates  higher  EWL  values  for  larger  charges  and  shallower  depths  within  this 
range.  Furthermore,  according  to  the  scaled  depths  listed  in  Table  1,  we  would  expect  similar  be¬ 
havior  from  the  case  W1D4  and  W2D5  (with  both  length  and  time  scaled  by  the  factor 

^2  « 1 .26 ).  This  is  clearly  reflected  in  Fig.  4,  as  well  as  the  similarity  between  the  W1D3  and 
W2D4  case.  This  scaling  is  not  exact  since  gravity  will  have  a  slightly  greater  influence  in  the 
larger  charge  case.  Plots  of  the  computed  values  EWL[t,Q,H)  for  H  =  A  ,6,  and  8  for  the  run  ma¬ 
trix  are  shown  in  Figure  5.  None  of  these  cases  produced  values  that  attained  the  desired  value  of 
4  ft  at  the  //  =  8  ft  height  level,  but  this  value  is  attained  at  //  =  4  ft,  for  charge  weights  IT  >  1 .5  . 


(U)  Plume  Velocity 

(U)  The  plume  velocity  may  also  enhance  its  protective  ability  as  a  barrier,  potentially  deflecting 
a  missile  by  imparting  a  vertical  momentum.  The  vertical  velocity  component  along  the  axis, 
v{t,Q,H) ,  is  shown  in  Fig.  6  for  the  same  cases  as  shown  for  the  EWL  values  in  Fig.  5.  Note  the 
large  negative  downward  velocities  reported  at  the  lower  //  =  4  ft  height  at  early  times.  These 
negative  velocities  indicate  the  downward  jet  seen  during  the  beginning  of  the  first  bubble  con¬ 
traction  shown  in  Fig.  2  at  times  t  =  0.1  and  t  =  0.2  seconds.  The  upward  moving  jet  can  be  ex¬ 
pected  to  attain  speeds  of  100  ft/s  or  larger,  during  the  time  of  the  first  bubble  expansion,  after 
which  it  is  decelerated  primarily  due  to  gravity.  We  remark  that  here  we  are  only  showing  the  ver¬ 
tical  velocities  along  the  axis,  and  are  thereby  ignoring  the  velocity  of  the  radial  plumed  ejected 
during  the  beginning  of  the  second  bubble  expansion. 
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Fig.  3.(U)  EWL  Contours  with  Plume  events  for  the  W2D3  case. 
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Fig.  4.  (U)  Comparison  of  EWL  contour  levels  for  the  original  run  matrix. 
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Fig.  5.  (U)  EWL  values  as  a  function  of  time  at  different  heights  for  the  original  run  matrix 
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Fig.  6.(U)  Vertical  plume  velocities  at  different  heights  for  the  original  run  matrix 
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(U)  Boundary  Effects 

(U)  The  effects  of  two  different  boundary  regions  were  also  compared  to  the  z  >  =  -20 ,  cases 

described  above.  For  these  comparisons  the  charge  weight  was  fixed  at  W  =  2  lbs,  and  the  initial 
depth  was  set  at  D  =  3  ft.  In  the  first  case,  the  initial  water  region  was  restricted  to  a  cylindrical 
“pool”  having  dimensions  r<R  =  9  and  0  s  z  s:  =  -8 .  The  purpose  of  this  simulation  was  for  a 
preliminary  investigation  for  designing  test  experiments  in  a  controlled  environment.  In  the  se¬ 
cond  case  only  the  bottom  boundary  was  changed  so  that  0  s  z  s:  =  -3 .  In  this  case  the  initial 
bubble  was  the  top  half  of  a  hemisphere  situated  at  the  bottom  boundary.  The  hemisphere  has  the 
same  volume  as  the  initial  bubble  for  a  spherical  charge  region.  In  each  case  the  grid  was  chosen 
to  match  the  original  fine  grid  with  uniform  cells  of  size  /i  =  0.05  for  z  ^  10  and  r  <  8 . 

(U)  Fig.  7,  exhibits  the  early  time  plume  and  bubble  contours  for  the  original  domain  compared  to 
the  two  restricted  domain  cases.  A  comparison  of  the  corresponding  EWL  contours  shown  in 
Fig.  8,  indicate  that  both  restricted  domains  increase  the  amount  of  water  thrown  upward.  In  the 
original  domain  the  water  surrounding  the  charge  is  pushed  outward  nearly  spherically  in  each 
direction.  In  the  case  of  shallow  water  when  the  charge  is  sitting  at  the  bottom  boundary,  all  of  the 
water  is  pushed  upward  or  radially  outward  (hemispherically).  Thus,  more  water  can  be  expected 
to  be  pushed  upward  than  in  the  original  domain  (since  the  total  volume  displaced  will  be  about 
the  same,  but  none  will  be  displaced  downward).  In  the  case  of  the  confined  pool,  a  significant 
portion  of  the  radial  velocity  is  converted  to  vertical  velocity  at  the  pool  side  boundaries,  again 
causing  more  water  to  be  directed  upward.  Based  on  these  observations,  we  conjecture  that  the 
case  of  a  shallow  plume  in  free  water  represents  the  case  of  lowest  EWL  values. 

(U)  SUMMARY  AND  CONCLUSIONS 

(U)  This  report  presents  a  numerical  investigation  of  the  amount  of  water  ejected  from  shallow 
depth  explosion  plumes  for  the  potential  use  as  a  barrier  from  incoming  projectiles.  The  computa¬ 
tional  code  BUB2D  was  used  for  the  simulations  based  on  the  generalized  hydrodynamics  model 
which  has  been  verified  and  validated  on  a  wide  variety  of  applications,  including  shallow  plume 
simulations  (see,  e.g.  [1,3-5]).  The  parameters  used  for  the  code  in  this  report  were  selected  such 
that  the  transition  from  a  non-venting  to  venting  bubble  occurred  at  the  approximate  scaled  depth 
value  C  =  DI  «  0.25 ,  consistent  with  values  reported  in  [8].  Using  these  same  parameters  and 
the  same  initial  conditions  as  the  shallow  plume  simulations,  a  grid  convergence  study  was  per¬ 
formed  on  the  spherically  symmetric  bubble,  demonstrating  the  errors  were  under  1%  for  the 
bubble  maximum  radius,  and  less  than  4%  for  the  first  period. 
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Fig.  13.  (U)  Early  time  boundary  effects  W2D3  case. 
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Fig.  14.  (U)  EWL  contours  for  the  W2D3  case  with  different  boundaries. 
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(U)  The  simulations  for  the  original  run  matrix  (eharge  weights  between  1  and  2  pounds,  and 
charge  depts.  Between  3  and  5  feet)  revealed  just  over  a  peak  equivalent  length  of  between  1  and 
3  feet  of  water  through  the  plume  center  is  ejected  when  measured  at  a  height  of  8  feet  above  the 
surface.  These  barriers  increased  to  between  2  and  5  feet  when  measured  at  a  height  of  6  feet,  and 
between  3  and  6  feet  at  a  height  of  4  feet  above  the  surface.  The  peak  plume  velocities  were  on 
the  order  of  100  ft/s,  but  these  decayed  after  the  time  of  the  first  bubble  maximum  volume  at  ap¬ 
proximately  t  =  0.1  seconds. 

(U)  The  effects  of  shallow  depths,  and  constrained  regions  were  shown  to  slightly  enhance  the 
EWL  levels.  In  the  first  case  this  demonstrates  that  if  an  effective  barrier  is  formed  in  deep  water, 
it  will  remain  effective  in  shallow  water,  provided  the  depth  is  not  so  shallow  that  the  bubble 
vents.  The  case  of  the  restricted  “pool”  shows  some  consideration  must  be  given  if  any  valida¬ 
tions  are  performed  in  a  confined  experimental  region. 


(U)  Future  Studies 

(U)  The  results  presented  here  represent  a  starting  point  for  the  feasibility  of  designing  an  effec¬ 
tive  protection  system  for  riverine  vessels.  There  are  several  critical  issues  and  studies  that  have 
not  yet  been  addressed  or  performed.  These  include 

1.  Further  model  validations.  The  best  validations  are  from  comparisons  to  actual  experi¬ 
ments.  However,  determining  EWL  values  accurately  is  not  trivial  as  demonstrated  in  the 
early  attempts  discussed  in  [11-13].  In  particular,  consistency  of  the  data  using  different 
measurement  techniques  was  not  adequately  demonstrated.  Alternatively,  the  current 
model  can  be  compared  to  other  models  which  have  also  been  validated  for  underwater 
explosion  problems.  One  such  code  is  DYSMAS  [14]  which  is  a  more  general  hydrocode 
than  BUB2D,  and  includes  a  compressible  equation  of  state  for  water,  thereby  treating 
shock  effects  directly. 

2.  Off  center  plumes.  The  EWL  values  reported  here  represented  the  best  case  since  they 
were  determined  by  integrating  in  a  straight  line  through  the  plume  center  (5  =  0 ).  The 
sensitivity  of  the  measurements  on  the  case  when  the  plume  is  not  in  the  direct  path  of  the 
plume  center  5  >  0  is  not  addressed  here. 

3.  Plume  velocity  effects.  The  velocity  measurements  were  done  for  the  purpose  of  future  in¬ 
cluding  the  momentum  imparted  on  the  incoming  missiles  due  to  the  vertical  velocity  of 
the  plumes.  Some  of  these  velocities  were  included  in  Fig.  6,  but  the  analysis  of  its  effect 
on  vehicle  protection  have  not  been  addressed. 

4.  Water  entry  effects.  The  water  barrier  concept  in  consideration  involves  the  entry  of  a 
charge  fired  into  the  water  between  the  vessel  and  incoming  missile.  The  water  entry 
event  is  expected  to  form  a  trailing  cavity  behind  the  projectile  which  will  probably  not 
close  before  the  charge  detonates.  The  effects  of  this  cavity  on  the  subsequent  bubble  and 
plume  formation  has  not  been  investigated. 

(U)  Addressing  these  topics  will  enable  the  reliable  assessment  of  the  water  barrier  protection  sys¬ 
tem.  With  accurate  simulation  tools,  system  performance  can  be  optimized,  without  the  need  of 
conducting  a  large  number  of  field  trials. 
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